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ABSTRACT 

The  generally  difficult  problem  of  a  system  with  time 
delay  is  studied  by  using  parameter  space  methods  and  com- 
puter analysis.   The  solution  is  proved  to  be  very  easily 
obtained  and  the  results  presented  are  extendable  under 
certain  circumstances  to  other  systems  and  useful  in  other 
fields  of  analysis. 
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I.   INTRODUCTION 

Dead  time  or  transport  lag,  defined  as  a  fixed  time  de- 
lay in  the  transmission  of  information  between  two  segments 
of  a  system,  is  a  phenomenon  that  occurs  frequently  in  every 
field  of  engineering.   The  synthesis  of  a  control  system 
where  dead  time  tends  to  occur  has  most  discouraging  aspects 
and  the  first  thoughts  of  the  engineer  confronted  with  the 
problem  should  be  how  to  go  about  changing  the  system  design 
so  as  to  eliminate  its  sources. 

However,  many  times  this  is  not  possible  due  to  physical 
limitations  in  the  system.   The  dead  time  may  be  due  to  hav- 
ing the  measuring  element  and  the  controller  actuator  neces- 
sarily separated  by  a  certain  distance,  as  in  the  extreme 
case  of  the  remote  control  of  a  space  vehicle.   In  other 
cases,  there  is  a. delay  between  the  time  at  which  the  con- 
troller actuator  is  moved  and  the  time  at  which  the  plant 
reacts  to  the  change,  as  in  the  case  of  an  engine-governor 
system. 

Dead  time  deteriorates  seriously  the  stability  of  a 
system  since  it  creates  a  phase  lag  that  constantly  increases 
with  co. 

The  analytic  treatment  of  the  problem  is  strenuous  and 
obtaining  an  exact  solution  is  very  difficult,  if  not  impos- 
sible.  The  Characteristic  Equation  (C.E.)  of  a  system  with 
time  delay  has  transcendental  nature  and  a  Maclaurin  series 
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of  the  exponential  term  shows  that  the  system  has  an  infinite 
number  of  roots. 

Bode  and  Root  Locus  analysis  are  practical  and  efficient 
methods  in  the  study  of  linear  systems,  but  in  the  case  of  a 
system  with  transport  lag  they  become  too  cumbersome.   They 
are  also  single  parameter  techniques.   However,  it  is  pre- 
cisely in  the  case  of  systems  with  transport  lag  that  the 
"Parameter  Space  Method"  (1),  (2),  appears  as  a  most  useful 
and  practical  aid.   This  method  has  application  whenever  the 
system  has  several  adjustable  parameters  and  consists  gener- 
ally in  the  mapping  of  the  S  plane  in  a  N  dimensional  space 
of  the  parameters  of  the  system.   With  this  mapping  completed 
it  is  possible  to  choose  the  N  parameter  values  that  allow 
the  system  to  have  a  certain  set  of  root  locations. 

As  previously  mentioned,  a  system  with  transport  lag  has 
a  C.E.  with  an  infinite  number  of  roots;  therefore,  if  we 
were  attempting  to  look  at  stability,  the  knowledge  of  the 
locations  of  some  of  the  roots  would  not  be  the  answer. 
Methods  are  available  for  testing  the  stability  of  systems 
with  transport  lag  (2)  ,  (3)  .   In  this  work  we  will  concen- 
trate on  finding  the  parameter  values  that  guarantee  a  near- 
optimal   behavior  of  the  system  recognizing,  however,  that 
analysis  of  stability  should  be  performed  after  this  synthe- 
sis phase  is  finished. 

The  study  of  the  general  problem  was  followed  by  appli- 
cation to  a  practical  example. 


II.   THEORETICAL  BACKGROUND 

Consider  the  two  block  diagrams  of  Figure  1.   The  trans- 
fer functions  of  the  two  systems  represented  are: 

,.,.,*  sififcCl! (1) 

X    1+G(S) .H(S) .e"ib 

T.F.2  -  ^ —  (2) 

■     1+G(S) .H(S) .e  ib 

We  see  that  the  form  of  the  C.E.  of  the  system  does  not 
depend  on  the  location  of  the  dead  time  block,  that  is,  the 
C.E.  is  the  same  independently  of  having  the  dead  time  block 
in  the  forward  or  backward  path.   This  allows  the  extension 
of  the  following  study  to  a  large  class  of  problems  where 
dead  time  occurs. 

The  term  exp(-TS)  gives  to  this  equation  a  transcendental 

nature.   In  fact,  the  Maclaurin  expansion  of  exp  (-TS)  is: 

— TS                  2       "3 
e     =  1  -  ST  +  ST   -  ST   +  (3) 

We  see  then  that  a  system  with  transport  lag  has  an  in- 
finite number  of  roots  which,  other  than  making  complete  sta- 
bility analysis  very  difficult,  has  the  effect  of  making  some 
graphical  methods,  as  Root  Locus  and  Bode  Diagram,  very  in- 
appropriate . 

The  stability  of  systems  where  dead  time  appears  is  al- 
ways a  factor  to  consider  since  stability  is  seriously 

-TS 
affected  by  it.   Consider  that  if  we  examine  e     through 
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sinusoidal  analysis  and  let  s  =  ju),  the  transport  lag  term 
becomes 

=  111  e-^  .    (4) 

We  see  then  that  the  transport  term  has  a  magnitude  of 
1  and  a  lag  phase  angle  of  coT .   The  phase  lag  added  by  the 
transport  lag  term  increases  as  co  increases.   The  system 
may,  therefore,  be  very  difficult  to  stabilize. 

However,  we  assumed  initially  that  the  nature  of  the 
system  is  such  that  all  the  roots  lie  in  the  left  half  of 
the  S  plane.   The  system  is  then  intrinsically  stable  and 
the  problem  is  to  find  the  values  of  the  adjustable  para- 
meters that  make  the  system  behave  in  an  optimal  manner  in 
relation  to  an  established  criterion.   We  concentrated  on 
finding  the  parameter  values  that  guarantee  the  existence 
of  a  set  of  dominant  roots  for  the  system,  that  is,  a  set 
of  roots  close  to  the  origin.   Stability  analysis,  using  for 
instance,  the  methods  in  references  (2)  and  (3),  should  be 
used  after  the  parameters  for  the  system  were  selected.   To 
determine  the  parameter  values,  the  parameter  space  method 
can  be  used  with  great  advantage  since,  as  will  be  shown, 
it  allows  us  to  relate  precisely  parameter  values  with  root 
locations.   We  considered  the  case  of  systems  with  three 
adjustable  parameters. 
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The  general  form  of  the  C.E.  of  a  system  with  transport  lag 
Is: 


N  N 

F(S)  «£  K.S1  +]T  K.SJ 

j=0 


-TS    n 
e     =  0 


(5) 


For 


i  =  0 
S  -  -O    +    jto 

-TS     Ta+jTw     10     ,  „   -        _  s   f,s 

e     =  e    J    =  e    (cos  Too  +  j  sxn  Tu)   (6) 

This  transformation  allows  the  separation  of  the  C.E. 
into  two  equations,  one  containing  the  real  and  the  other 
the  imaginary  part.   That  is: 

F(a,oj)  =  F  (a, to)  +  jF  (a,oo)  =  0 


implies  that : 


F  (a,oj)  =  0 
e 

F  (a,aj)  =.0 
o   ' 


(7) 
(8) 


If  we  consider  that  the  system  has  three  adjustable  parameters 
A,  B  and  C,  then  they  appear  in  F  (a,oj),  and/or  F  (a, to). 
Fixing  one  of  the  parameters  (for  instance,  parameter  A), 
equations  (7)  and  (8),  may  be  solved  for  the  two  parameters 
B  and  C.   Recapitulating,  we  have  that  equations  (7)  and  (8) 
are  functions  of  the  variable  parameters  A,  B  and  C  of  O    and 
to  from  the  assumed  roots  -O   +  jto.   Fixing  all  these  variables 
but  two  of  the  parameters,  for  instance,  B  and  C,  simulta- 
neous solution  of  the  two  equations  allows  the  determination 

of  the  two  values  B   and  C  ,  that  together  with  the  fixed 

o       o  ° 

value  of  A   guarantees  that  the  system  has  a  pair  of  conju- 
gate roots  at  -a   +  jw  . 

o  —    o 
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Repeating  this  process  for  all  different  values  of  oo  in 
a  certain  range  permits  the  determination  of  different  set- 
pairs,  B  C  ,  as  in  Figure  2.   Of  course,  there  is  an  infinite 
r     *   o  o 

number  of  these  solutions.   The  union  of  these  solutions 

B  C   is  a  line  in   the  parameter  plane  of  axis  B  and  C  as 
o  o 

is  shown  in  the  Figure  3  (a  line).   Working  similarly,  but 

with  the  variation  of  a  instead  of  oj,  it  is  possible  to 

obtain  similar  lines  in  the  parameter  plane  B-C  for  which  co 

is  fixed  (w  lines) . 

The  intersection  of  a  certain  line  of  fixed  a)  (to   line) 

o   o 

and  another  line  of  fixed  a   (a   line)  gives  a  point  in  the 

o    o        ° 

parameter  plane  from  which  it  is  possible  to  read  the  pair 

it        * 

of  parameters  B   C    that  guarantee  the  existence  of  a  pair 

o   o 

of  roots  for  the  system  of  the  form  -O      +  id)  .   See  Figure  4. 

'  o  —  J  o  b 

Repeating  the  above  process  for  different  values  of  A 

allows  the  determination  of  a  line  (or  a  group  of  points) 

in  the  three-dimensional  space  with  ordinates  A,  B  and  C  for 

which  every  point  corresponds  to  assigning  the  system  the 

pair  of  roots  -a   +  i  co  .   See  Figure  5. 

o  —  J  o  6 

If  we  assume  that  the  system  has  a  third  dominant  root 
-6  -  of  a  real  nature,  we  can,  fixing  the  value  of  A,  deter- 
mine a  line  of  fixed  6(6       line)  in  the  parameter  plane  B-C 
as  in  Figure  6.   This  is  so  because  if  we  substitute  the 

solution  6   in  the  characteristic  equation  and  fix  its  value 
o  n 

together  with  the  value  of  A,  then  this  equation  becomes  the 
equation  of  a  line  in  the  plane  B-C. 
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Assuming  different  values  for  A,  but  keeping  6  constant, 
the  repetition  of  this  process  for  each  value  of  A  results 
in  a  group  of  lines  in  which  <5  is  constant,  and  these  lines, 
if  this  repetition  were  performed  for  an  infinite  number  of 
values  of  A,  form  a  surface  (surface  6  )  in  the  parameter 
space  A,  B  and  C  for  which  the  system  has  a  real  root  6  , 
as  in  Figure  7 . 

In  the  three-dimensional  parameter  space,  the  pierce 

point(s)  of  the  line  -a   +  iw   with  the  plane  -6   gives  one 
r  o  —    o  r        o° 

or  more  points  that  correspond  to  the  set(s)  of  parameters 

A  ,  B   and  C  ,  for  which  the  system  has  a  set  of  three  roots 
o*   o       o* 

-0   +  iw   and  6   that  we  assumed  would  be  the  dominant  ones, 
o  —  J  o       o 

Notice  that,  since  we  are  mostly  interested  in  finding 
roots  that  lie  close  to  the  origin,  the  range  of  variation 
of  a,  a),  and  6  during  the  process  of  mapping  the  S-plane  in 
the  parameter  space  does  not  need  to  be  very  large. 

Of  course,  such  an  extensive  manipulation  of  equations 
and  plot  of  curves  requires  the  use  of  a  computer. 

We  will  proceed  now  to  apply  these  theoretical  considera- 
tions to  the  study  of  a  practical  problem,  specifically,  to 
the  determination  of  the  set  of  near-optimal  parameters  to  be 
used  in  the  operation  of  an  engine-governor  system.   This 
problem  was  suggested  by  the  reading  of  reference  (4)  and 
subsequent  exchange  of  correspondence  with  one  of  its  authors 

In  this  paper,  it  is  suggested  that  the  determination  of 
the  parameter  settings  must  be  made  by  assuming  that  the 
dominant  roots  of  the  system  occupy  a  certain  well-defined 
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position,  specifically  that  the  real  and  imaginary  parts  of 
the  complex  roots  must  have  equal  magnitude  equal  also  to 
the  magnitude  of  the  third  (real)  root. 

The  parameter  space  method  was  used  to  find  new  combina- 
tions of  parameter  values  that  would  permit  us  to  determine 
if  the  root  locations,  as  assumed  in  the  aforementioned 
paper,  would  correspond  to  an  optimal  solution,  and  to 
search  for  new  solutions.   In  addition,  the  results  permit 
study  of  related  problems  such  as  the  problem  of  sensitivity. 


III.   PARAMETER  PLANE  STUDY  OF  AN  ENGINE-GOVERNOR  SYSTEM 

A.   GENERAL  CONSIDERATIONS 

Most  gas  engines  exhibit  a  certain  degree  of  dead  time, 
that  is,  an  interval  of  time  between  a  change  in  fuel  flow, 
and  the  corresponding  change  in  driving  torque.   This  time 
lag  is  made  of  a  number  of  components  of  which  the  most 
important  are:  the  delay  in  "sensing"  the  speed  change  that 
follows  the  load  change;  the  delay  within  the  governing 
mechanism  during  which  the  reaction  of  the  sensing  element 
is  being  translated  into  governor  action;  the  time  required 
for  the  governor  action  to  correct  the  fuel  injection  pump 
setting;  the  delay  that  must  take  place  before  the  corrected 
fuel  change  can  be  introduced  into  an  engine  cylinder;  the 
time  required  to  convert  the  corrected  fuel  change  into 
engine  torque. 

Dead  time  is  found  to  be  detrimental  to  stable  operation 
and  tends  to  increase  speed  transient  errors.   It  also  deter- 
mines the  maximum  rate  at  which  the  engine-governor  combina- 
tion can  be  operated.   The  purpose  of  this  study  is  to  find 
the  best  isochronous  governor  parameters  so  that  optimum 
transients,  for  the  operation  of  gas  engines  of  several 
types  be  obtained. 

The  study  of  ideal  governor  settings  in  the  presence  of 
dead  time  will  be  made  in  the  parameter  space.   This  method 
leads  to  a  computer  process  for  analyzing  the  variable  system 
parameters.  The  process  is  restricted  to  a  two  parameter  space 
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B.  ENGINE  ASSUMPTIONS 

The  driving  torque  incremental  change  AQ  is  related  to 
the  small  variations  of  the  governor  actuator  AZ  by  the 
equation : 

AQ  =  C6  e"TS  AZ  (9) 

-TS 
where  e     is  the  dead  time  operator. 

The  minimum  amount  of  dead  time  T  that  can  be  expected 
is  given  by  the  following  formulae: 

T  =  20/N   for  a  Diesel  engine  (10) 

T  =  40/N   for  a  two-cycle  gas  engine  (G.E.)(11) 
T  =  60/N   for  a  four-cycle  gas  engine (G . E. ) (12) 
in  which  N  is  the  rotational  speed  in  RPM . 

The  assumptions  made  relative  to  the  engine  rotating 
parts  are  that  the  accelerating  torque  is  equal  to  the  dif- 
ference between  the  load  torque  and  the  driving  torque  as 
expressed  by  the  following  equation: 

~  S  AN  -  AQ  -  AL  (13) 

Combination  of  (9)  and  (13)  gives: 

^  S  AN  =  C,  e~TS  AZ  -  AL  (14) 

Change  to  a  system  of  per  unit  variables  gives  the 

engine  equation  used  in  the  analysis: 

T   Sn  =  e"TSZ  -  I  (15) 

M 

C.  GOVERNOR  ASSUMPTIONS 

It  is  assumed  that  the  isochronous  pressure  compensated 
type  of  governor  is  used.   The  differential  equation 
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describing  its  dynamics  is: 

(S2  +  a1S)Z  =  -V2(S  +  a3)n  (16) 

where  u2  =  (Nr/Zr)a2. 

D.  ENGINE  GOVERNOR  COMBINATION 

Combination  of  equations  (15)  and  (16)  gives: 

(S3+  a,S2+  ~  e~TSS  +  -4"1  e"TS)n  =  -  jjr-S(S+a-)£    (17) 
M  M  M 

E.  USE  OF  THE  PARAMETER  SPACE  METHOD 

We  start  with  the  engine-governor  equation  (Equation  17) 

(S^+a.sN-  -Ae    ibS  +  -f-^-e  ib)n  =  -  i_S(S+a.H 
M  M  M 

From  equation  (17),  the  C.E.  can  be  extracted: 

•5        2    ^2   -TS      ^2a3   -TS 
SJ  +  a  SZ  +  ^   e  ibS  +  -f-±   e  ib  =  0  (18) 

M  M 

Since  we  have  three  adjustable  parameters  in  the  charac- 
teristic equation,  we  can  control  the  position  of  three  of 
its  roots.   We  assumed  then  that  the  system  had  three  domi- 
nant roots  similarly  to  what  would  happen  if  the  delay  in 
the  same  equation  were  equal  to  zero.   For  these  three  domi- 
nant roots  we  had  two  choices  with  respect  to  their  intended 
location,  that  is,  we  could  choose  to  have  three  real  roots, 
or  two  complex  conjugate  roots  and  a  real  one.   Because  we 
were  interested  in  having  a  small  rise  time  for  the  tran- 
sients of  the  system,  we  chose  the  second  case,  that  is,  we 
assumed  that  the  system  should  exhibit  a  pair  of  complex 
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conjugate  roots  and  a  single  real  root  close  to  the  origin 
in  the  S  plane. 

In  accordance  with  this  choice  we  made: 

S  =  -c  +  ju) 

We  have  then  that: 


2      2     2 
S   =»  (a   -w)-j  (2au)) 


S3   =    (3aco2    -   a3)    +  j(3a2oo   -   w3) 
-TS 


The  exponential  term  e     becomes: 


-TS     -T(-a+jco) 
e     =  e       J 

TO       -jwT 
=  e    e  J 

Ta 


=  e    (cos  wT  -  j  sin  u)T) 
Substituting  (19),  (20)  and  (21)  in  (18)  gives: 


(19) 
(20) 


(21) 


2    3        2      3        2     2 

(3aco  -  a  )  +  j  (3a  a)  -  a)  )+  a, (a   -  to  -  j2aoo) 


M2   Ta 
+  —  e   (cos  coT-j  sin  u)T)  (-a+j  w) 

M 

y2a3   Ta 
+  — = e   (cos  U)T-j  sin  wT)  =  0 

M 


(22) 


Both  the  real  and  the  imaginary  parts  of  this  equation 
must  be  equal  to  zero: 


2  ~3 


2  ,.2 


2  Ta 
Saco^-a^+a,  (az-a)z)+  — e   (a>  sin  coT-a  cosa)T+a„cosuT)=  0   (23) 

M  J 


3a  a)-w3-2a1aw+  =— e   (a  sinuT+u  cosuT-a„sina)T)  =  0       (24) 

M 

Equations  (23)  and  (24)  form  a  system  that  can  give  solu- 
tions for  any  two  parameters  if  the  others  are  fixed.   If  we 
take  ex..  and  y9  for  parameters  of  the  plane,  and  fix  the  value 
of  a^,  we  can  get  a  complete  representation  of  the  loci  of 
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the  points  fa  (w)  ,  y_  (u>)  ,  a  =  constant  or  0^(0),  y2(cr), 
(0  =  constant!  for  which  the  C.E.  has  a  solution  of  the  form 
-a  +  ja). 

The  curves  determined  as  above-curves  a  and  U)   are  loci 
of  all  parameter  pairs  (a,  and  y„)  for  which  the  solutions 
-0*  +  ja)  and  -a-jid  exist. 

If  we  denote  by  6  the  dominant  real  root  of  the  C.E.,  we 
can  by  the  same  methods  get  a  plot  of  6  lines  in  the  a,  ,  y« 
parameter  plane. 

We  can  repeat  the  process  of  finding  a,  to  and  6  curves 
for  as  many  values  of  a_  as  we  want. 

For  a  certain  a.,  =  constant,  the  intersection  of  selected 

a   and  w   curves  is  the  loci  of  values  of  the  parameters  a, 
00  1 

and  y_  for  which  the  selected  solution  -a   +  iw   exists. 
2  o  —  J  o 

In  the  three  dimensional  orthogonal  space  (a..  ,  y?  and 

a„)  the  loci  of  pairs  a,  and  y~  is  a  line  along  which  it  is 

guaranteed  that  any  readings  of  a, ,  y«  and  a„  guarantee  the 

existence  of  the  above  form  of  solution.   In  a  similar  way, 

the  set  of  lines  of  constant  6  ,  one  for  each  value  of  a„ 

o  3 

forms  a  plane  in  the  three  dimensional  space.  «• 

The  single  or  multiple  intersections  of  the  aforementioned 

line  with  the  plane  of  constant  6=6   gives  the  set(s)  of 

parameters  a,,  y„  and  a„  for  which  the  system  has  a  real 

root  6   and  a  pair  of  conjugate  roots  -o   +  iw  .   However, 
o        r  J  b  0—0 

since  we  are  primarily  interested  in  finding  the  near-optimal 
set  of  parameters  for  the  governor,  that  is,  the  set  of  para- 
meters that  offers  a  "better"  transient,  we  will  not  be 
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looking  for  a  set  of  the  three  parameter  space  parameters  as 
indicated  above,   Instead,  we  will  fix  the  location  of  the 
pair  of  complex  roots,  which,  as  explained  above,  will  allow 
the  determination  of  a  line  in  the  three  dimensional  para- 
meter space.   We  will  then  proceed  to  test  the  system  with 
parameters  corresponding  to  points  of  this  line.   Of  course, 
the  parameters  corresponding  to  each  of  these  points  deter- 
mine a  different  location  for  the  dominant  real  root. 

A  computer  program  was  written  to  plot  the  a,  a)  and  <5 
curves  on  a  parameter  plane  with  axis  a.,  and  y  ?  for  a  number 
of  values  of  the  parameter  a„  (computer  program  #1).   This 
program  was  designed  with  the  purpose  of  making  it  adaptable 
to  studies  of  any  type  of  engine-governor  combination  that 
may  be  described  by  equation  (17). 

We  performed  the  computer  analysis  first  for  a  hypothe- 
tical system  in  which  the  transport  lag  was  zero  and  later 
for  each  of  the  three  practical  cases  considered — Diesel 
engine,  two-cycle  gas  engine,  and  four-cycle  gas  engine- 
each  of  them  with  the  same  C.E.,  but  with  different  dead 
times  as  in  equations  (10)  ,  (11)  and  (12)  .   The  problem  of 
obtaining  parameter  values  that  guarantee  certain  root  loca- 
tions is  thus  solved.   We  still  needed,  however,  any  kind  of 
way  of  testing  the  systems  with  each  of  the  sets  of  para- 
meters obtained  through  the  parameter  space  computer  program. 
We,  naturally,  made  use  of  a  second  computer  program  (compu- 
ter program  //2)  which  was  designed  under  the  following  heading. 
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IV.   TESTING  OF  THE  ENGINE-GOVERNOR  SYSTEM 

We  want  a  program  to  test  each  of  the  systems  with  dif- 
ferent groups  of  parameter  settings.  From  equation  (17)  we 
determine  the  transfer  function  (T.F.): 

-  |  S(S  +  a1) 

f  -       " (25) 

J,    c2^  y2  -TSC    U2a3  -TS 
S  +  a..S  +  — e    S  +  — e 

M  M 

The  denominator  of  the  above  T.F.  is  rearranged  as 
follows : 


3       2    y2   -TSC    U2a3   -TS 
S   +  a,S   +  —  e    S+  — e 

M  M 

2  ^2   -TS 

=  S  (S  +  a  )  +  —  e     (S  +  a  ) 

x     M  J 

-TS 
y    e    (S  +  a.) 

=  SZ(S  +  a,) (1  + 


'1  '  v      T      2 
1       XM   SZ(S  +  c^) 


Then,  the  T.F.  becomes: 

1 


S(S  +  ax) 

n  _  M 

y    e    (S  +  CO 

SZ  (S  +oO  (1  +  ~  -5 — 

x  XM   S   (S  +  a  ) 

M  G 


-TS,„  .     v   1  +  GH 

y9   e 
1  + 


y2 

e          (S    +    a3) 

T 
M 

S2(S    +    a    ). 

G    = 

1 

T    S 
M 

(26) 


This  implies  that:      G  =  -  ~~  (27) 


ye    (S+a  ) 
S  (,S+a, ; 
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Using  block  diagram  description,  we  can  partition  G  and 
H  in  the  following  form: 


G  = 


H  = 


l/a3s+l 
1/a-S+l 


^2a3 
a. 


And  the  complete  block  diagram  of  the  feedback  system  becomes 
as  represented  in  Figure  8.   This  configuration  allows  the 
use  of  C.S.M.P. 
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V.   PROCESSING  OF  A  PRACTICAL  PROBLEM 
AND  ANALYSIS  OF  RESULTS 

The  design  of  governors  requires  the  choice  of  a  "stan- 
dard" type  of  response  to  load  changes.   Since  a  change  in 
load  unavoidably  causes  a  transient  change  in  speed,  the 
usual  application  has  the  following  requirements: 

a)  Recovery  should  be  fast  and  with  negligible  differ- 
ence in  steady  state  speed; 

b)  The  maximum  deviation  from  steady  state  speed  should 
be  less  than  some  specified  percent; 

c)  The  transient  recovery  should  not  overshoot  the 
steady  state  value. 

Trade-offs  between  these  requirements  may  be  necessary  in 
some  applications. 

For  a  system  with  no  time  delay,  but  with  third  order 
characteristics,  the  best  solution  for  items  a)  and  c)  is 
a  set  of  two  complex  and  one  real  root,  such  that  the  real 
root  is  numerically  equal  to  the  real  parts  of  the  complex 
roots,  and  such  that  the  imaginary  and  real  parts  of  the 
complex  roots  are  equal.   For  typical  applications,  a  numeri- 
cal choice  which  gives  a  satisfactory  solution  is  a  real 
root  at  -3  and  complex  roots  at  -3  +  j3.   (A  study  of  re- 
sults based  on  other  choices  is  made  later.)   Justification 
of  this  is  given  by  Figures  9,  10,  11  and  12,  which  show  the 
transient  responses  for  a  system  with  no  time  delay,  with 
complex  roots  at  -3  +  j3  and  with  several  locations  of  the 
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real  root.   Notice  that  Figure  11  is  the  one  that  corresponds 
to  having  the  real  root  6  equal  to  -3  and  that  Figures  9  and 
10  don't  satisfy  item  a)  while  Figure  12  shows  an  overshoot 
which  goes  against  the  requirement  posed  in  item  c) . 

We  wanted  to  find  the  near-optimal  parameter  settings 
for  the  governor  described  in  reference  4  for  each  of  the 
types  considered  in  the  same  paper.   We  assume  in  all  cases 
that  the  steady  state  speed  of  the  engine  is  450  RFM. 

We  used  some  practical  data  as  it  appears  in  the  same 
reference : 

TM  =  5.88 

T   =  KR/N  =  20/450  for  a  Diesel  engine 
=  40/450  for  a  2  cycle  G.E. 
-  60/450  for  a  4  cycle  G.E. 

We  needed  also  to  know  what  was  the  behavior  of  a  system 
governed  by  equation  (17),  but  with  zero  delay.   This  could 
be  helpful  in  establishing  the  effect  of  the  time  delay  and 
was  used  in  determining  parameter  values  that,  entered  as 
data  in  computer  program  #2,  allowed  the  plotting  of  Figures 
9,  10,  11  and  12. 

In  each  case  the  parameter  plane  program  was  run  for 
values  of  the  parameter  a,,  equal  to  0.5,  0.6,  0.7,  0.8,  0.9, 
1.0,  1.1,  1.2,  1.3,  1.4,  1.5,  1.55,  1.6  and  1.7. 

We  decided  in  all  cases  that  this  system  should  have  a 
pair  of  roots  at  -3  +  j3  and  so  we  concentrated  on  finding 
the  location  of  the  point  of  the  parameter  planes  where  the 
curves  to  =  3.0  and  a  =  3.0  intersected.   Had  the  choice  of 
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root  locations  been  different,  we  could  still  have  used  the 
same  parameter  plane  plots  since  they  contain  curves  for 
values  of  a  and  a)  equal  to  0.5,  1.0,  1.5,  2.0,  2.5,  3.0,  3.5, 
4.0,  4.5,  5.0,  5.5  and  6.0.   The  reason  why  we  did  not  run 
the  program  for  values  of  a„  higher  than  1.7  was  that  we 
found  that  in  each  case  for  values  of  a_  higher  than  a  well- 
determined  value,  the  O    =  3.0  and  oj  =  3.0  curves  on  which  we 
concentrated,  would  not  cross  in  the  parameter  plane.   This 
is  justified  by  the  following: 

For  example,  in  the  case  of  the  four  cycle  gas  engine 
with  a  =  w  =  3.0,  we  have  from  equation  (23)  that: 

0  4 
54  +  y2  x  0.17  x  e    (3  sin  0.4-3  cos  0.4  +  a.cos  0.4)=  0 


.254  y2  (0.92  a3  -  1.59)=-54 


-212.6 


2   -92a3  -  1.59 


Assuming  that  the  physical  nature  of  the  parameters  re- 
quires that  they  have  positive  values  has  the  effect  that 
y„  must  be  positive,  which  implies  that 

„  1.59  ~  .  _0 
a3  <  —92  "  1'73 

The  analysis  of  the  parameter  plane  plots  proceeded  as 
follows:  for  each  a~  =  const.,  we  determined  the  intersection 
of  the  a  =  3.0  and  w  =  3.0  curves  and  read  the  corresponding 
coordinate  values  of  a,  and  y?. 

Plotting  the  point  correspondent  to  these  parameter 
values  in  page  2  of  the  parameter  plane  plot,  (page  of  curves 
of  constant  6)  allowed  the  determination  of  the  location  of 
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the  real  root  for  the  considered  value  of  a».   This  value  of 
6  was  further  checked  by  using  computer  program  //3. 

Tables  1,  2,  3  and  4  list  in  each  row  the  parameter 
values  obtained  by  this  process,  and  that  correspondingly 
must  be  used  in  testing  the  dynamics  of  each  of  the  systems 
so  that  they  present  a  pair  of  complex  roots  at  -3  +  j3, 
and  the  value  of  the  real  root  determined  as  described. 

These  tables  were  used  to  plot  Figures  13,  14  and  15 
that  represent  the  projections  in  each  of  the  orthogonal 
planes  that  form  the  parameter  space  of  the  points  of  this 
space  that  guarantee  that  each  of  the  systems  will  have  a 
pair  of  roots  at  -3  +  j3.   Notice  in  these  plots  that  the 
points  correspondent  to  a  certain  system  are  disposed  "in  a 
line".   Notice  also  that  due  to  this  form  of  disposition, 
if  we  were  to  work  with  a  system  described  by  equation  (17), 
but  with  a  delay  different  from  those  considered  here,  the 
determination  of  parameter  values  that  would  give  the  same 
root  locations  would  be  possible  by  simply  using  interpola- 
tion. 

Figures  16  through  55  are  some  of  the  parameter  plane 
plots  obtained  by  using  the  parameter  plane  program  in  the 
IBM  360/65  digital  computer.   These  are  only  a  part  of  the 
large  number  of  parameter  plane  plots  that  we  had  to  obtain 
to  construct  Tables  1,  2,  3  and  4.   We  chose  those  that 
either  give  parameter  values  close  to  those  found  to  the 
"best"  or  the  ones  that  give  real  "bad"  parameter  sets. 


30 


a3 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.49 

1.5 

1.55 

1.6 

1.7 


al 
6.6 
6.7 
6.9 
7.1 
7.2 
7.4 
7.5 
7.9 
8.3 
8.5 
8.97 
9.0 
9.2 
9.3 
9.8 


126.9 
132.3 
138.6 
145.8 
150.0 
158.0 
167.0 
177.0 
188.0 
198.2 
209.5 
210.6 
216.9 
225.0 
242.0 


•0.60 
■0.75 
•0.91 
■1.09 
■1.29 
■1.49 
•1.68 
■1.95 
•2.28 
•2.54 
•3.00 
3.03 
3.26 
•3.33 
•3.82 


a3 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 
1.1 
1.2 
1.3 
1.4 
1.42 
1.5 
1.55 
1.6 
1.7 


al 

7-6 

7.8 

8.2 

8.5 

8.7 

9.0 

9.5 

10.0 

10.3 

10.9 

11.05 

11.7 

12.0 

12.4 

13.5 


142.2 
144.5 
148.5 
154.8 
165.0 
175.0 
187.0 
201.6 
217.0 
235.0 
238.6 
253.8 
264.0 
279.0 
306.0 


■0.60 
■0.77 
•0.95 
•1.15 
■1.29 
■1.58 
•1.88 
■2.18 
■2.43 
■2.82 
■3.00 
•3.45 
■3.68 
•3.85 
•4.88 


Table    1    -    T    =    0 


Table    2    -    T    =    20/450 
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a3 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.49 

1.5 

1.55 

1.6 

1.7 


al 
6.6 
6.7 
6.9 
7.1 
7.2 
7.4 
7.5 
7.9 
8.3 
8.5 
8.97 
9.0 
9.2 
9.3 
9.8 


126.9 
132.3 
138.6 
145.8 
150.0 
158.0 
167.0 
177.0 
188.0 
198.2 
209.5 
210.6 
216.9 
225.0 
242.0 


-0.60 
-0.75 
-0.91 
-1.09 
-1.29 
•1.49 
•1.68 
■1.95 
■2.28 
■2.54 
■3.00 
•3.03 
3.26 
■3.33 
•3.82 


a3 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 
1.1 
1.2 
1.3 
1.4 
1.42 
1.5 
1.55 
1.6 
1.7 


al 

7.6 

7.8 

8.2 

8.5 

8.7 

9.0 

9.5 

10.0 

10.3 

10.9 

11.05 

11.7 

12.0 

12.4 

13.5 


142.2 
144.5 
148.5 
154.8 
165.0 
175.0 
187.0 
201.6 
217.0 
235.0 
238.6 
253.8 
264.0 
279.0 
306.0 


•0.60 
■0.77 
■0.95 
■1.15 
■1.29 
•1.58 
•1.88 
•2.18 
2.43 
■2.82 
3.00 
3.45 
3.68 
3.85 
4.88 


Table    1    -    T    =    0 


Table    2    -    T    =    20/450 
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a3 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 

1.34 

1.4 

1.5 

1.55 

1.6 

1.7 


al 
9.5 
9.7 
10.2 
10.7 
11.3 
11.9 
12.8 
13.8 
14.6 
15.26 
16.1 
17.8 
18.9 
20.3 
23.8 


150.3 
157.5 
167.0 
182.0 
196.0 
213.3 
234.0 
256.5 
286.0 
301.4 
320.4 
367.2 
390.0 
430.2 
514.8 


-0.63 
-0.78 
-1.02 
-1.19 
-1.44 
-1.70 
-2.04 
-2.50 
-2.72 
-3.00 
-3.36 
-3.87 
-4.25 
-4.72 
-6.80 


a3 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.55 

1.6 

1.7 


al 

13.4 
12.7 
13.4 
15.4 
16.8 
18.4 
21.0 
24.1 
28.5 
36.8 
48.4 
65.3 
89.5 
421.0 


193.0 

207.0 

225.0 

251.0 

280.0 

317.0 

368.0 

438.0 

532.0 

707  .0 

971.0 

1300.0 

1859.0 

9092.0 


-0.67 
-0.90 
-1.02 
-1.25 
-1.58 
-1.91 
-2.15 
-2.54 
-3.00 
-3.69 
-4.50 
-4.95 
-5.50 
-7.20 


Table    3    -    T    =    40/450 


Table    4    -    T    =    60/450 
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Each  of  the  3  systems  was  simulated  including  their  res- 
pective time  delays  and  the  parameter  values  as  in  each  row 
of  tables  1-4.   Figures  56  through  70  are  plottings  of  the 
transients  of  each  of  the  systems  when  subjected  to  100%  of 
load  change.   We  obtained  the  transient  responses  for  every 
set  of  parameters  in  those  tables,  but  we  are  including  only 
as  many  as  needed  to  illustrate  the  following  analysis. 

Figure  71  is  a  plot  of  the  transients  of  each  of  the 
systems  for  the  case  in  which  the  dominant  real  root  has  a 
value  of  -3.0.   We  can  see  that  all  these  transient  res- 
ponses are  very  similar.   Thus  the  three  roots  considered 
in  the  cases  of  the  systems  in  which  the  delay  is  non-zero, 
are  really  dominant  over  all  the  others,  and  the  determina- 
tion of  the  near-optimal  parameters  for  a  system  with  trans- 
port lag  consists  of  finding  the  parameter  values  that 
determine  the  location  of  the  dominant  roots  to  be  the 
ideal  location  of  the  roots  for  a  system  with  the  same 
equation,  but  with  zero  time  delay.   (This  criteria  is 
restricted  to  a  certain  portion  of  the  S  plane  as  will  be 
seen  later . ) 

Examination  of  Figure  71  shows  that  for  each  of  the 
systems,  the  settling  time  is  very  approximately  the  same; 
the  differences  in  time  delay  determine  only  the  amplitude 
of  the  transient,  this  being  proportional  to  the  amplitude 
of  the  delay,  which  makes  sense  since  the  greater  the  delay, 
the  more  time  the  governor  takes  to  interpret  the  change 
and  set  the  needed  corrections. 
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A  comparative  analysis  of  Figures  56  through  70  shows 
that  the  transients  of  each  system  with  parameter  settings 
that  correspond  to  having  the  real  root  closer  to  the  origin, 
have  larger  settling  times,  while  in  the  opposite  case,  the 
recoveries  are  faster,  but  an  overshoot  appears. 

The  study  of  the  method  of  determination  of  parameter 
sets  is  completed.   We  tried  to  make  use  of  the  large  amount 
of  data  available  by  proceeding  with  the  analysis  of  some 
other  subjects.   The  first  of  these  relates  to  the  choice 
of  location  of  the  pair  of  complex  roots. 

For  a  linear  system,  we  know  that  the  same  root  pattern, 
that  is,  having  the  real  root  numerically  equal  to  the  real 
part  of  the  complex  roots  and  this  real  part  equal  in  magni- 
tude to  the  imaginary  part,  should  give  the  same  type  of 
transient,  except  of  course  for  time  scaling.   In  other 
words,  the  three  groups  of  dominant  roots, 

r  )  •    c )  -  C"  ) 

\-3+j3/  \-4+j4/  \-10+jl0/ 

should  all  give  identical  transients  for  a  linear  system, 
except  for  the  time  scale.   However,  in  the  case  of  a  system 
with  time  delay,  this  should  not  apply,  and  there  should  be 
some  numerical  limit  on  the  pattern,  depending  on  the  value 
of  T .   In  other  words,  if  we  tried  to  make  this  root  pattern 
fall  far  to  the  left  of  the  S  plane,  either  the  parameters 
would  become  negative-therefore  physically  unrealizable-or 
some  other  roots  would  become  dominant,  or  would  even  fall 
in  the  right  half  of  the  S  plane.   The  answer  to  this  should 
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be  looked  at,  because,  given  that  this  pattern  were  realiz- 
able more  to  the  left  of  the  S  plane,  certainly  the  tran- 
sient would  become  faster,  and  the  magnitude  of  the  maximum 
deviation  could  also  possibly  be  smaller. 

We  started  by  finding  the  sets  of  parameters  that  would 
give  to  the  Diesel  governor  system  root  sets  respectively 
at  (-4  and  -4  +  j4),  and  (-5  and  -5  +  j5).   The  determina- 
tion of  parameter  sets  went  as  before,  that  is,  in  each  case, 
we  looked  in  page  one  of  each  of  the  plots  of  const,  a   for 
the  intersections  of  the  a  =  4  with  a)  =  4  and  0    =    5    with 
0)  =  5 .   For  each  of  these  intersections,  we  determined  in 
page  two  of  the  plot  the  corresponding  value  of  the  real 
root  6,  and  were  finally  able  to  determine  the  values  of 
the  parameters  a,,  a„  and  y„,  that  guarantee  the  desired 
root  locations.   The  parameter  sets  determined  are  displayed 
in  the  following  table: 


DOMINANT 
ROOTS 

al 

a3 

^2 

-4,  -4+j4 

16.16 

1.84 

444.5 

-5,  -5+j5 

20.91 

2.25 

643.0 

We  used  computer  program  //2  to  determine  the  responses 
of  the  system  while  working  with  these  sets  of  parameters, 
and  we  obtained  the  transients  displayed  in  Figures  72  and 
73.   Looking  at  Figure  72,  it  can  be  seen  that,  with  the 
parameters  set  such  that  the  system  has  dominant  roots  at 
-4  and  -4  +  j4,  the  speed  transient  is  faster  and  of  smaller 
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amplitude  than  any  one  of  those  previously  found,  which 
guarantees  that  the  correspondent  choice  of  parameters  is 
closer  to  the  "optimal".   However,  the  response,  displayed 
in  Figure  73  of  the  system  with  dominant  roots  at  -5  and 
-5  +  j5,  is  not  good  under  item  c)  of  the  established  cri- 
teria.  Here  the  recovery  is  very  fast  and  the  amplitude 
of  the  transient  is  even  smaller  but  an  overshoot  shows. 
This  proves  that  our  hypothesis  about  the  existence  of  a 
limit  in  the  process  of  bringing  the  root  pattern  more  to 
the  left  of  the  S  plane,  was  reasonable. 

Another  topic  of  interest  was  the  one  related  to  toler- 
ances in  the  parameter  values  and/or  sensitivity  of  the 
systems  to  parameter  change.   The  possibility  of  and  ease 
in  performing  such  studies  is  one  of  the  great  advantages 
of  the  parameter  space  method,  since  the  parameter  plane 
plots  display  sets  of  curves  that  allow  to  relate  changes 
in  parameter  values  with  changes  in  root  location.   Figures 
74  and  75  show  a  parameter  plane  plot  where  we  marked  the 
point  that  confers  to  the  engine-governor  system  a  set  of 
roots  at  -A  and  -A  +  jA.   The  -.-  lines  correspond  to  varia- 
tions of  a  certain  percent  in  the  determined  values  of  the 
parameters  a.,  and  y~  .   The  possible  changes  in  root  locations, 
due  to  changes  in  parameter  values,  is  quite  obvious.   For 
example,  the  considered  variation  in  the  value  of  a. ,  while 
^2  is  kept  constant,  can  cause  variation  of  03  between  D  and 
E  and  variations  in  the  values  of  a  between  B  and  C. 
Obviously,  similar  considerations  apply  to  page  2  of  each 
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plot  and  allow  the  study  of  changes  in  the  value  of  6.   The 
analysis  of  the  effects  of  changes  in  a~  is  a  little  more 
difficult,  since  it  requires  the  use  of  plots  correspondent 
to  adjacent  values  of  a  .   Here  we  should  plot  the  point 
correspondent  to  the  nominal  values  of  a1  and  y~  in  the 
adjacent  plots.   The  point  plotted  in  the  adjacent  plots 
shows  immediately  the  changes  in  a,  w  and  6  due  to  the  varia 
tion  of  a_. 

Of  course,  if  it  were  desirable  to  do  so,  the  whole 
process  could  be  reversed;  that  is,  we  could  easily  find 
out  the  new  values  of  the  parameters,  if  it  were  required 
to  change  the  location  of  the  roots  in  relation  to  the 
values  set  initially. 

Finally,  notice  that  the  parameter  plane  plots  are  a 
great  aid  in  studies  of  sensitivity  of  the  systems  since 
the  slope  of  the  lines  plotted  gives  direct  measure  of  that 
characteristic. 
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VI.   CONCLUSIONS 

The  association  of  the  parameter  plane  method,  with 
appropriate  computer  analysis,  converted  a  rather  difficult 
problem  into  an  almost  elementary  one.   In  effect,  after 
passing  the  difficulties  in  writing  and  debugging  the  para- 
meter plane  program,  the  analysis  of  the  information  avail- 
able is  so  easy,  that  an  ill-advised  observer  may  not  notice 
all  the  trouble  involved  in  the  solution  of  the  same  problem, 
if  it  were  attempted  by  other  methods. 

The  problem  of  finding  the  set  of  near  optimal  parameters 
with  time  delay  consists  of  finding  the  parameter  values 
that  guarantee  that  the  dominant  roots  of  the  C.E.  of  the 
system  will  be  located  according  to  a  special  pattern.   The 
quality  of  the  response  improves  with  the  displacement  of 
the  dominant  roots  to  the  left  of  the  S  plane  but,  this  pro- 
cess has  a  limit  depending  on  the  value  of  the  delay.   The 
determination  of  this  limit  at  which  the  quality  of  the  res- 
ponse should  be  near-optimal  has  to  be  made  through  a  trial 
and  error  process. 

The  method  described  in  this  study  is  very  easy  to  follow, 
and  the  information  available  allows  for  investigation  in 
many  areas  of  interest,  such  as  sensitivity  of  the  systems  to 
variations  in  parameter  values,  analysis  of  permitted  toler- 
ances in  these  values,  optimization  of  performance  under  a 
defined  cost  function,  and  others. 
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Figure  1.   Control  systems  with  transport  lag  in  the 
forward  and  backward  paths. 
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3.    3 


Figure  2.   Plot  of  simultaneous  solution  of  equations 
(7)  and  (8). 
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Figure  3.   O       line  in  the  A  =  const,  plane. 


Figure  A. 


a   and  oo   lines  in  A  =  const,  plane. 
o       o  r 
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Figure  5.   a  ,  u   line  in  A,  B,  C  parameter  space 
o    o 
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Figure  6.   a   line  in  A  =  const,  plane, 

b  o 


Figure  7.   6   plane  in  A,  B,  C  parameter  space 
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Figure  9.   Transient  response  of  a  hypothetical  system 
with  zero  time  delay. 
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Figure  10.   Transient  response  of  a  hypothetical  system 
with  zero  time  delay. 
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Plot  of  6  lines. 
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Plot  of  6  lines. 
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Plot    of    6    lines. 
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Plot    of    a,     u)    lines. 
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Plot    of    a,    a)    lines. 
T=  2o/4fO         A3=  /2 
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Plot  of  6  lines . 
T=  20/450  A3=  1.2 
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2o/4fO        A3=  1J 
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Y   SCALE=90.0   UNITS/INCH 


Fig.  20 
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Plot    of    6    lines. 
T=  20/450         A3=  /J 
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Y  SCALE=90.0  UNITS/INCH 
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X  SCALE=06.0  UNITS/INCH 
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X  SCALE=06.0  UNITS/INCH 
Y  SCALE=90.0  UNITS/INCH 
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Plot    of    a,     to    lines." 

T=2o/4S0         A3=/5" 

X   SCALE=06.0    UNITS/INCH 
Y   SCALE=90.0   UNITS/INCH 


Fig. 32 
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Plot  of  6  lines . 

T=  20/4SO        A3=  IS 
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Fig. 3? 


oi: 


"wr 


67 


Plot    of    a,     to   lines. 
T=  20/45-0         A3=  f,7 

X    SCALE=06.0   UNITS/INCH 
Y   SCALE=90.0   UNITS/INCH 


Fig.  34 
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Plot    of    6    lines. 
T=20/4fO        A3=/7 

X  SCALE=06.0  UNITS/INCH 
Y  SCALE=90.0  UNITS/INCH 


Fig.  35" 
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Plot    of    a,     03    lines. 
T=4o/45"o  A3  =0.5" 

X   SCALE=06.0   UNITS/INCH 
Y   SCALE=90.0   UNITS/INCH 


Fig.  3G 
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Plot    of     6    lines. 
T=  40/450  A7>=0S 

X   SCALE=06.0   UNITS/INCH 
Y   SCALE=90.0   UNITS/INCH 


Fig. 37 


OJO 


ex, 


71 


X   SCALE=06.0   UNITS/INCH 
Y   SCALE=90.0   UNITS/INCH 


Fig. 36 
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Plot    of     6    lines. 
T=  40/45-O        A3=/3 

X  SCALE=06.0  UNITS/INCH 
Y  SCALE=90.0  UNITS/INCH 


Fig.  39 
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Plot    of    6    lin< 
T=40/4f0        AZ=/.4 
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T=40/ '450       A3=/7 
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Y   SCALE=90.0   UNITS/INCH 


Fig.  42 
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Plot    of    6    lines. 

i=4o/4$o       A3=  17 

X  SCALE=06.0  UNITS/INCH 
Y  SCALE=90.0  UNITS/INCH 
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Plot    of    a,    a>    lim 
T=6o/4fO  AZ=0.'f 
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Fig.  44 
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Plot  of  6  lines. 
T^GO/45'O        A  3=  0.5" 

X  SCALE=06.0  UNITS/INCH 
Y  SCALE=90.0  UNITS/INCH 
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Plot    of    a,    a)    lim 
1=  &0 1450         A3=/2 
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Plot  of  6  lines. 
T=&V '4S0  A3=/2 

X  SCALE=06.0  UNITS/INCH 
Y  SCALE=90.0  UNITS/INCH 


Fig.  47 
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X  SCALE=06.0  UNITS/INCH 
Y  SCALE=90.0  UNITS/INCH 
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Plot    of    a,     03    lines. 
T=C,0/ffO         A3=/4 

X   SCALE=06.0   UNITS/INCH 
Y   SCALE=90.0   UNITS/INCH 
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Plot    of    6    lines. 
1=&0/45'0        A3=/4 

X   SCALE=06.0   UNITS/INCH 
Y   SCALE=90.0   UNITS/INCH 

Fig.  57 
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Plot    of    a,     a)    lines. 

T-60/4S&        A3- 7.5" 

X  SCALE=06.0  UNITS/INCH 
Y  SCALE=90.0  UNITS/INCH 


Fig. 52 
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Plot    of    6    lines. 

1=Go/4SO       kl=/S 

X   SCALE=06.0   UNITS/INCH 
Y    SCALE=90.0   UNITS/INCH 
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Plot  of   a,  tn  lines. 

T=6o/4fO         A3=/7 

X   SCALE=06.0   UNITS/INCH 
Y   SCALE=90.0   UNITS/INCH 
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Plot    of    6    lines. 
T=  &0/4S0         A3=/-7 

X   SCALE=06.0   UNITS/INCH 
Y   SCALE=90.0   UNITS/INCH 
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Figure  71.   Speed  transient  of  the  four  systems  when  the 
three  dominant  roots  are  located  at  -3  and  -3  +  j3. 
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Fig.  74 
Roots  displacement  versus 
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